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ABSTRACT 


The thermal effects in the Concentric Canister “missile” Launcher (CCL) due to a 
fire in an adjacent compartment are simulated using computational fluid dynamics (CFD). 
A commercial code developed by CFD Research Corporation (CFDRC) was used to 
implement the process. This study developed a model for the center section of the A- 
module, placed in the forward missile launcher of the Arleigh Burke (DDG-51) class 
destroyer. Two fire scenerios are applied to the aft bulkhead of the launcher. The first 
fire scenerio 1s indicative of the high temperature fire caused by unburned Exocet missile 
propellant experienced by the USS STARK (FFG-31) in the Arabian Gulf. The second 
fire scenerio applied to the model simulates the conditions caused by a ruptured shipboard 
F-76 diesel fuel tank, due to a collision at sea. The fire scenerios are used to predict the 
time and location of the critical cook-off temperatures of the missile’s solid and liquid fuel 


propellants in the CCL. 
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I. INTRODUCTION 


A. DESCRIPTION OF THE PROBLEM 

In the Navy’s quest to improve the survivability of the Vertical Launch System 
(VLS) on U.S. combatants, Naval Surface Warfare Center Dahlgren Division (NSWCDD) 
has developed the Concentric Canister Launcher (CCL) as an alternative to the current 
Mk 41 system. Ref 1] 

The CCL is a ready-to-fire self-contained system that can launch diverse munition 
types and sizes at any angle (Figure 1). Enhanced survivability measures in the design 
include a shock collar and passive armor protection. The shock collar is mounted in close 
proximity to the main deck to allow effective dissipation of an underwater blast away from 
the keel of the ship. The passive armor protection is provided by two titanium concentric 
canisters that are connected by interior anti-fragmentation shields. This combination of 
shields and strong inner and outer shells provides greater survivability of the missiles from 
enemy weapon penetration into the launcher. 

NSWCDD has test fired the NATO Sea Sparrow, Standard Missile (SM-2 Blk IV) 
and a navalized-version of the Army’s Advanced Tactical Missile (ATACM) from the 
CCL. They have also conducted extensive finite element analyses of imposed stresses on 
the steel shock collar and CFD simulations of the heat and mass transfer of the missile on 
the hemispherical lower cap of the canister. [Ref. 2] 

To date, NSWCDD has not tested the effects on the CCL magazine due to high 
temperatures caused by a fire in an adjacent compartment. The first step in fire modeling 


the CCL system was thesis research performed by Callaham (1996). He conducted a one- 
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Figure 1. AA-Size Concentric Canister Weapons Module. “From Ref. [4].” 
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dimensional lumped parameter analysis on a single canister located next to the adjacent 
hot boundary [Ref. 3]. The fire modeling of the entire launcher is an important step to 


ensure the CCL meets present survivability criteria. 


B. BACKGROUND 

Currently, the Mk 41 VLS installed on the Arleigh Burke (DDG-51) class 
destroyer represents the state of the art in vertical launch design. The DDG-51 has two 
VLS that are designed to hold Tomahawk (T-LAM B, C, D), Standard (SM-2 Blk I, III, 
IV) and the Vertical Launch Anti-submarine Rocket (VLA). The forward magazine can 
carry up to 32 missiles and the aft magazine can carry up to 64. Table 1 shows that the 
drawback of carrying a large number of missiles onboard is that each contains several 
hundred pounds of energetic material (i.e., solid propellants, liquid fuels and explosives). 
A fire in the VLS involving this much energetic material could be catastrophic and result 


in the loss of the ship. 
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Table 1. Munitions Contained in VLS Loadout. “From Ref. [4].” 







In 1987, the USS STARK (FFG-331) was struck by two Iraqi Exocet missiles 
while on patrol in the Arabian Gulf. The missiles failed to detonate upon impact, 
however, the missiles’ solid propellant continued to burn at very high temperatures. 
Figure 2 shows that although the fire was contained in the midships section, the forward 
Mk 13 Standard Missile magazine still sustained thermal damage. 

The STARK had state of the art fire protection systems, yet due to the high 
temperatures and damage to parts of these systems, they proved inadequate in fighting a 
“weapon-induced” fire [Ref. 4]. This type of fire is the most dangerous to the Burke class 
destroyers. The unburned Exocet fuel produced fires with temperatures in excess of 
3000° F (1922 K). The minimum temperature required to melt the steel used in Navy 
ships is roughly 2800° F (1811 K). 

In 1992, the USS INGERSOLL (DD-990) was involved in a major collision at sea 
with a Singapore merchant vessel. Four of INGERSOLL’s forward fuel oil storage tanks 
were ruptured, which immediately flooded two of the forward compartments with fuel oil. 
One of the spaces flooded was the five inch gun magazine, while the other was a space 
immediately adjacent to this compartment. Fortunately, a fire in INGERSOLL was 
avoided. The crew maintained a thick layer of firefighting foam over the top of the two 
compartments and shored much of the damaged starboard side, allowing the ship to 


maintain ten knots enroute to the Singapore shipyard. 
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Figure 2. USS STARK (FFG-31). Starboard Side View Summary of Damage. 
“From Ref, [5].” 





C. PREVIOUS WORK 


Bowman and Lee developed a one-dimensional heat transfer model to predict the 
time-temperature profiles in the missile canister of a Mk 41 VLS for the DDG-51. Their 
scenerios used a fire in an adjacent compartment, which increased temperatures from 
ambient conditions to 2000°F (1366 K) in a period of five minutes, and then maintained 
that temperature for the duration of the fire (Figures 3 and 4). [Ref 4] This fire scenerio 
is indicative of the F-76 fuel oil fire that INGERSOLL potentially faced during her 
collision in the Straits of Malacca. 

The goal of the analysis was to determine the approximate time for the missile’s 
propellants to reach their critical cook-off temperatures. The cook-off temperature for the 
solid propellants is 300° F (422 K) and 460° F (511 K) is the ignition temperature for the 
liquid fuel. Bowman and Lees’ model predicted that the missile motors and warheads 
were in danger of self-ignition in approximately 50 minutes, if more than one adacent 
space was on fire. For a fire in a single adjacent space, the canister temperature was 
predicted to reach only 245° F (377 K) after 900 minutes. Although this temperature is 
unlikely to cook-off the solid propellants and warheads, the Otto II fuel in the VLA may 
be at risk. [Ref 4] 

Callaham generated a MATLAB 4.2.c computer algorithm that predicted the time 
temperature profile in a single CCL canister containing a Tomahawk missile. Using four 
different fire scenerios, he varied the temperature and heat fluxes applied to the adjacent 
bulkhead of the CCL compartment to determine the critical value of the cook-off 


temperature of the propellant at the geometric center of the missile. [Ref 3] 
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Figure 3. Ship’s Structure at VLS Boundary. “From Ref. [4].” 
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Figure 4. Canister Time Temperature Profile for Various Fire Sizes. 
“From Ref. [4].” 


Callaham’s model accounted for a single canister located in the center of the 
nearest row to the heat affected adjacent bulkhead. He also performed a simulation on a 
canister located at the corner of the nearest row. Figure 5S shows the thermal resistance 
network of the CCL canister. The model did not, however, account for the reflected 
radiation received from two of the four bulkheads, the upper and lower steel supports and 
canisters in the compartment, as well as, the circulating natural convection on the side of 
the canister that faced away from the hot boundary. 

The most severe scenerio analyzed was one modeled as a fire fueled by a missile’s 
residual solid propellant, as experienced by the STARK. The maximum temperature 
imposed was 3000° F (1922 K) for a rate of increase of 97.3° F/sec (54 K/sec). Table 2 
shows the resultant times for the center of the canister to reach the solid propellant and 
liquid fuel cook-off temperatures for the STARK scenerio. Table 2 also shows the 
resultant times for the scenerio used by Bowman and Lee, applied by Callaham to his 
computer algorithm. Based on his data, he concluded that under the worst case fire 
scenerio of the STARK, the ship’s personnel have approximately 45 minutes to gain 


control of the fire or cool the weapons in the magazine prior to cook-off. [Ref. 3] 
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Table 2. Cook-Off and Ignition Times for Callaham Model. “After Ref. [3].” 
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Figure 5. Thermal Resistance Network of the CCL Magazine. “From Ref. [3].” 


D. OBJECTIVES 


The purpose of this study is to expand upon the previous one-dimensional lumped 
parameter analysis performed on vertical launch missile magazines, adjacent to a 
compartment affected by a high temperature shipboard fire. The new three dimensional 
model includes all of the magazine’s bulkheads and the addition of several canisters inside. 

The initial step in developing such a model was to generate a meshing code and 
map it to a finite difference algorithm. The mesh generation code needed to account for 
the rectangular geometry of the launcher enclosure, as well as, the cylindrical geometry of 
the canisters and the accompanying radiation shape factors. To develop an accurate finite 
difference algorithm, the effects of surface radiation, natural convection, and conduction 
needed to be included in the program. The construction of each of these codes alone was 
considered to be overly complex given the existence of commercial CFD codes. 

The decision was made to procure a commercial CFD program that included a 
mesh generation pre-processor, as well as, a powerful three dimensional algonthm for 
solving the continuity, momentum and energy equations. The baseline model is 
representative of the midsection of an A-module CCL compartment, installed in the 
forward missile magazine of the DDG-51. Two separate fire scenerios are applied to the 
affected boundary condition, indicative of a high temperature fire caused by unburned 
Exocet missile propellant or shipboard F-76 diesel fuel. The data generated by the 
program can provide a visual representation of the time-temperature profile within the 
CCL compartment, as a result of a fire in an adjacent compartment, to aid in the 


development of optimum firefighting systems and procedures for the launcher. 


Il. COMPUTATIONAL FLUID DYNAMICS PROGRAM 


A. OVERVIEW 


The CFD program used to analyze the thermal effects in the CCL was developed 
by CFD Research Corporation (CFDRC). The program incorporates three separate, yet 
interactive, codes titled CFD-GEOM, CFD-ACE and CFD-VIEW (Figure 6). 

CFD-GEOM is the pre-processor geometric modeling and mesh generation code 
that enables the user to develop the initial architecture of the problem. CFD-ACE allows 
the user to select the phenomena to be included (1.e., heat transfer, turbulence, radiation), 
the values for the material properties and initial conditions, and the type of boundary 
conditions used (i.e., adiabatic, isothermal, inlet etc.). This code also allows the user 
limited control of the algonthm by enabling the choice of the number of iterations, the type 
of differencing scheme and the amount of constraint, or relaxation, that the solution can 
change as it 1s iterated through each cycle of the nonlinear governing equations. The final 
code, CFD-VIEW, is a post-processor data visualization tool that enables the user to 


picture the output through colored planes, vectors or animation. 


B. THEORY 


CFD-ACE uses a finite volume approach to solve the Favre-averaged Navier 
Stokes (FANS), continuity and energy equations for velocity, temperature, pressure and 
density. The geometry is divided into a number of finite volumes, or cells, each generated 


by a structured grid that uses single I, J, K indices to identify a particular gridpoint and 
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Figure 6. CFDRC Computational Fluid Dynamics Computer Code Overview. 
“From Ref. [6].” 
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direction within the finite volume. Each set of cells comprises a domain, which can be 
combined into a set of sud-domains to facilitate an easier construction of the overall 
modeled geometry. Since the CCL 1s a simple parallelepiped with internal cylinders, a 
cartesian coordinate system is used. The following ts a brief summary of the discretization 
of the governing equations and numerical iterative process used to solve for the unknown 
variables, as addressed in the CFDRC Theory Manual. A more in-depth discussion of the 
theory can be found in the manual and its reference section. [Ref. 6, 7] 

In CFD programs, the partial differential equations (PDE) that govern the 
transport of flow quantities are solved using numerical techniques. The techniques involve 
the discretization of the PDEs on a computational grid, the formation of a set of algebraic 
equations, and the solution of these equations at each grid point through an iterative 
algorithm. 

The nature of the flow problem in the simulation of thermal effects in the CCL 
involves turbulent fluid flow and conduction, convection and radiation heat transfer. The 
standard fluid flow governing equations in tensor form are: 


The continuity equation 


O 
Ft Be (pu) =0 


where p is the fluid density and u; 1s the j" Cartesian component of the instantaneous 
velocity. 


The conservation of momentum equation 


oO Oo. _ _P , os 
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where p is the static pressure, f; is the body force and tj is the viscous stress tensor 
defined as 


Ou: 


wa = (Oe 4 Oey _ 2 (ae 
aren 


eam 
3 on 
where up is the fluid dynamic viscosity and 6; 1s the Kronecker delta. 


And the static enthalpy equation 
O O 
ee ae eu = ty t+ 


where h is the enthalpy and q; is the } component of the heat flux defined as 


OT 
ioe -K— 
q ae 


where K 1s the thermal conductivity. 
All of the governing equations possess a common form and can be generalized by 


the equation 


op 


me 


0 0 mror 
By Pe) GT a (Ir 


where the symbol @ may represent any of the velocity components, enthalpy or the 
turbulent kinetic energy and dissipation rate. S is the source term and I 1s the diffusion 
coefficient. For the energy equation, the diffusion coefficient is the effective 


diffusivity defined as 


rat ,w_R 
Oo Ort Cp 


where c, is the specific heat, o is the Prandtl number and the subscript t represents 


turbulent quantities. 


The governing equations are discretized using a finite volume approach. As 
mentioned above, the model is divided into discrete sub-domains, or cells. CFD-ACE 
employs a co-located cell-centered variable arrangement where all flow variables and fluid 
properties are stored at the cell center P. 

The following equation is derived for each variable in each finite volume 

Ade = Ache + Awhw + Andy + Ashs + angu + ard. +S 
where P is the center location and E, W, N, S, H and L represent values of the center of 
the six adjacent cells. The variable ‘a’ represents the effects of convection and diffusion 
across the cell faces. 

Using the east face (E) as an example, the convection and diffusion discretized 
terms are written as follows: 

Convection Term 
—peAcuede 
Diffusion Term 


_p, (oe) A; 
Ar 


The source term, S, contains terms other than the convection and diffusion effects, 
such as the pressure gradient, boundary conditions and other algorithm dependent effects 
(i.e., under-relaxation). Ag is the area across the east face and Ag is the distance from the 
cell center to the center of the adjacent face. 

The differencing scheme determines how the cell face values are calculated. The 


differencing scheme used in this model is the first-order upwind scheme. In this scheme, 
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de is taken to be the value of at the upstream grid point, which depends on the flow 


direction at the cell face E: 


de= op ifug 20 
Oe = OF if up < 0 


The continuity equation is not written in the form of the general convection- 
diffusion equation, as were the conservation of momentum and energy equations. Since, 
the pressure and velocity terms are strongly coupled in the momentum equation, the 
pressure gradient term is accounted for differently than the other source terms. 

For each iteration, the pressure and velocity terms will equal the current value plus 


some unknown correction according to the following equation: 


p=p*+p 
u=u*+u 


When the correction approaches zero then we can be reseasonably confident that 
the value of pressure and velocity are accurate for the given flow field. The values of 
these corrections are determined by using the momentum equation to develop a functional 
relationship between the pressure and velocity fields. The pressures are then substituted 
for the velocities into a discretized continuity equation, and the simultaneous matnx of 
equations is solved over the entire flow field. The solution procedure is summarized in 
Figure 7. NITERP is the number of iterations for the continuity equation. NTSTEP is the 
number of time steps involved. NITER is the number of iterations of the system of 
equations for each time step. 

Figure 7 shows that for each iteration, within a transient time step, a system of 
equations must be solved for each dependent variable. A benefit to using the structured 


grid geometry is that it produces a banded matrix of coefficients for each dependent 
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Figure 7. CFD-ACE Flow Solver Solution Procedure. “After Ref. [6].” 
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variable. The solution methods for each system of equations are Backward Euler and a 
version of the Forward Differencing Technique. 

In order to ensure that the solution for each iteration does not diverge, under- 
relaxation must be applied to the dependent and auxiliary variables. Under-relaxation 
constrains the amount that each variable can change from one iteration to the next. The 
dependent variables (u, v, w, k, ¢, H) are modified using an Inertial Factor. A linear under- 


relaxation technique 1s applied to the auxiliary variables (p, T, pL, p). 


II. MODEL 


A. GEOMETRY 


1. Selection of CCL Module 


Since the CCL Concept is under development as an alternative weapons system 
for the Surface Combatant of the 21° century (DD-21), the DDG-51 Flight II destroyer 
was used as the platform for construction of the CCL model in CFD-GEOM. The A- 
module standard size was selected since it is designed to hold 35 Tomahawk missiles and 
thus fit in the DDG-51’s forward missile launcher. Figure 8 shows the dimensions of the 
A-module and an equidistant distribution of 35 Tomahawk missiles inside. 

In order to develop a credible fire scenerio for the forward launcher, the fire was 
placed in a single compartment, aft of the forward launcher, to simulate the same 
conditions experienced by the STARK mass conflagration (Figure 9). The six 
compartments, adjacent to the common bulkhead of the CCL enclosure, are assumed to 
occupy an entire deck, without any intermediate longitudinal bulkheads or partitions to 


further contain the fire. 


2. Grid Distribution 


The simulations were carried out through the use of a Silicon Graphics Indigo XS- 
4000 workstation with 96 megabytes of RAM and an external 9 gigabyte hard drive to 
supplement the internal hard drive. The limitations of the computations using this 
workstation and the CFD-ACE version 2.0 program, updated as of 01 March 1997, were 


65,000 grid points and 500 boundary conditions per domain. 
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Figure 8. A-Module Standard Size with 35 Tomahawk Missiles. 
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Figure 9. Sideview of the DDG-51 Forward Missile Launcher. 
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The initial simulations on the A-module involved a geometric model of the entire 
launcher, including the correct dimensions and material specifications of the HY-80 
horizontal steel supports. The canister geometry initially included two concentric titanium 
cylinders, separated by air, and an internal aluminum core. This overall geometry 
exceeded the allowable number of grid points and boundary conditions, and was therefore 
reduced to the center one-third of the missile launcher, located between the two steel 
supports, with the canister geometry changed to include only an outer cylinder of solid 
titanium. The height of the center section is 2.33 meters. 

The geometry included in the model was reduced for the following reasons. Initial 
simulations on the entire missile launcher, involving only a single mode of heat transfer 
(i.e., conduction or convection) showed that the sections of the canisters located adjacent 
to the fire were directly affected by the thermal effects, as expected, and those sections 
that were separated by a steel support received only minor residual effects. Therefore, 
since the fire was located in the center compartment aft of the launcher, the decision was 
made to model only the center section of the CCL. Additionally, this selection was made 
because it tended to accurately model the location of the solid and liquid propellants of the 
missiles contained in the launcher. The selection of solid canisters was based entirely on 
the conservation of nodal points, and the determination that the effects between the 
canister walls were minimal in comparison to the large scale of the CCL model. 

Based on the updated reduced model of the CCL, the first simulations were 
conducted on a launcher that contained only the single aft row of canisters immediately 
adjacent to the affected fire boundary. Figure 10 shows the distribution of nodal points 


along the edges of this model. The vertical edges each contained 15 nodal points and the 


Ze 





transverse edges each contained 34. 

A second row of canisters was then added directly behind the first row to gain an 
understanding of how the buildup of canisters affected the row immediately adjacent to the 
affected fire boundary. Figure 11 shows the distribution of nodal points along the edges 
of this model, with the vertical and transverse edges containing the same number of nodal 
points as Figure 10. 

In the construction of the single- and double-row models, a high level of accuracy 
was not anticipated for the initial simulations due to the small number of nodal points, and 
hence large nodal point spacing. A crude method used to verify the accuracy of the 
output data was to vary the number of nodal points in the longitudinal direction. A single 
direction was chosen in order to maintain consistency throughout the computer runs. The 
longitudinal direction was the one chosen because it reflected the greatest change in the 
temperature gradients. 

Figure 12 shows the distribution of nodal points in the longitudinal direction for 
the single- and double-row models. Note, the grid points within the canisters were 


monitored by the CFD-ACE flow solver program and thus never altered. 


B. THERMOPHYSICAL MODEL 


Once CFD-GEOM has determined the grid and boundary condition locations and 
generated files containing this information, the files are then read into the CFD-ACE code. 
The CFD-ACE flow solver code allows the user to choose thermophysical properties, 
modes of heat transfer, boundary condition values, turbulent flow models and numerical 


solution methods. 
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Figure 12. The Total Number of Nodal Points on Each Longitudinal Edge. 
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The goal of the computations was to input enough data into the flow solver 
algorithm to accurately visualize the thermal effects within the CCL and to determine the 
elapsed time until propellant cook-off, while keeping the model simple, so as to limit the 
time required to run each simulation. Therefore, the canisters were modeled as solid 
titanium cylinders and the steel bulkheads and supports were given an emissivity 
equivalent to HY-80 steel [Ref. 8]. Table 3 lists the material property information used in 


all of the CFD-ACE simulations. 


DENSITY 4930 kg/m 
ar ee 
a 
| 
i 









_— 83 MOLECULAR WEIGHT 29 kg / kmol 


PRANDTL NUMBER 


—— 


Table 3. Thermophysical Property Data 
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Each simulation was carried out as a turbulent incompressible flow with surface 
radiation, conduction and non-Boussinesq natural convection heat transfer. Table 4 lists 
additional information used in the CFD-ACE simulations. The input file is contained in 


Appendix C. 
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Table 4. CFD-ACE Input Data 





C. FIRE SCENERIOS 


The simulated fire scenerios applied to the CCL model represent two of the most 
severe mass conflagrations that face our warships today. In order to realistically apply 
each of these boundary conditions, the maximum temperature and rate of temperature 
increase from ambient conditions must be included in the input file of the computer 


program. 
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The first scenerio analyzed is indicative of the severe inferno fueled by an anti-ship 
missile’s unburned solid propellant, as experienced by the STARK in 1987. The maximum 
temperature used to model this fire is 3000° F (1922 K), and climbs rapidly from ambient 
conditions to this maximum temperature in 30 seconds, where it is maintained for the 
duration of the simulation (Figure 13) [Ref. 3]. 

The second scenerio represents a less severe, though more likely to occur, fire 
fueled by a ruptured or leaking fuel oil storage tank, as faced by the INGERSOLL in 
1992. The maximum temperature used to model this fire 1s 2000° F (1366 K), and climbs 
from ambient to this maximum temperature in a relatively slower rate of five minutes [Ref. 
3]. Figure 14 shows the plot of this fire scenerio, in addition to the points along the curve 


that are added to the CFD-ACE input file. 
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Figure 13. Unburned Exocet Missile Rocket Propellant Fire Scenerio. 
“After Ref. [3].” 
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Figure 14. Common Shipboard Fuel Fire Scenerio. “After Ref. [3].” 
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IV. RESULTS 


This study models the thermal effects on the center section of a single-row and 
two-row CCL module using the CFD-ACE program. An unburned Exocet missile 
propellant fire and an F-76 shipboard fuel fire are used to simulate the scenerios along the 
affected boundary. The remaining three walls and upper and lower steel supports are 
modeled as adiabatic boundaries since the amount of heat conducted through them is 
negligible to the amount of heat transfer within the center section. 

Once the geometric, thermophysical and flow models are determined in the CFD- 
ACE and CFD-GEOM programs, the remaining challenge is to decide on the length of 
time for the processor to run and the size of the time step to use. Using the nodal 
distribution of round | (see Figure 12), a comparison was made between the following 
four runs: 

1. Time Interval: 0 - 60 seconds, Time Step: 1 second 

2. Time Interval: 0 - 600 seconds, Time Step: 10 seconds 

3. Time Interval: 0 - 1800 seconds, Time Step: 10 seconds 

4. Time Interval: 0 - 3600 seconds, Time Step: 60 seconds 

The decision was made to use a time interval of 600 seconds (10 minutes) with a 
time step of 10 seconds for all of the simulations. The larger time intervals and time step 
produced poor residual plots and ran for an extensive amount of time. None of the 
simulations reached the critical cook-off temperatures in less than one minute, therefore 


this time interval was determined to be too short. 


Although a time step of 10 seconds may have been judged to be too long for an 
accurate solution, the amount of time that these simulations took under this condition was 
roughly 3 to 4.5 days of continous processing. Decreasing the size of the time step would 
have been unproductive to this study. 

CFD-ACE allows the user to input a set of locations that the code monitors for 
changes in selected variables during each simulation (see Appendix C). The value is 
monitored at the center of each cell, and is identified in the input file by its domain and the 
indices of the top left nodal point of the cell. An example of the code used in the input file 
can be found at the end of Appendix C. 

In the CCL simulations, the temperatures within the center canister of the row 
adjacent to the fire boundary were monitored at a height of 0.69 meters (2 nodal points) 
above the bottom support, at three internal locations. This height was determined to 
roughly estimate the location of the missile’s propellant within the canister. Figure 15 
shows the resultant times to reach the critical cook-off temperatures of the solid and liquid 
propellants for rounds three and five of the Exocet and F-76 fire scenenos. The time- 
temperature profiles for all five rounds of the two fire scenerios can be found in 
Appendices D through M. 

Additional aspects of the temperature profiles at a height of 0.69 meters are also 
generated using the CFD-VIEW program. One profile uses a line probe that is sliced 
through the center canister of the front row, just left of centerline, to display the 
temperature as a function of longitudinal displacement from the fire boundary. Figures 16 
through 20 show these images for the Exocet fire scenerios and Figures 21 through 25 


show the F-76 fire scenerios. 
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Figure 15. Elapsed Times to Cook-Off for CCL Simulation Rounds Three and Five. 
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Figure 16. Temperature vs. Longitudinal Index, Displaced One Nodal Point Left of 
Centerline for Round 1 of the Exocet Fire Scenerio. 
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Figure 17. Temperature vs. Longitudinal Index, Displaced One Nodal Point Left of 
Centerline for Round 2 of the Exocet Fire Scenerio. 
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NOTE: Forward Skin of Missile is at x = 19. Aft Skin of Missile is at x = 22. 
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Figure 18. Temperature vs. Longitudinal Index, Displaced One Nodal Point Left of 
Centerline for Round 3 of the Exocet Fire Scenerio. 
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Figure 19. Temperature vs. Longitudinal Index, Displaced One Nodal Point Left of 
Centerline for Round 4 of the Exocet Fire Scenerio. 
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Figure 20. Temperature vs. Longitudinal Index, Displaced One Nodal Point Left of 
Centerline for Round 5 of the Exocet Fire Scenerio. 
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Figure 21. Temperature vs. Longitudinal Index, Displaced One Nodal Point Left of 
Centerline for Round 1 of the F-76 Fire Scenerio. 
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Temperature (K) 


NOTE: Forward Skin of Missile is atx = 11. Aft Skin of Missile is at x = 14. 
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Figure 22. Temperature vs. Longitudinal Index, Displaced One Nodal Point Left of 
Centerline for Round 2 of the F-76 Fire Scenerio. 
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Temperature (K) 


NOTE: Forward Skin of Missile is at x = 19. Aft Skin of Missile is at x = 22. 
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Figure 23. Temperature vs. Longitudinal Index, Displaced One Nodal Point Left of 
Centerline for Round 3 of the F-76 Fire Scenerio. 
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NOTE: First Row: Forward Skin of Missile is at x = 11. Aft Skin of Missile is at x = 14. 





Second Row: Forward Skin: of Missile isatx= 5. Aft Skin of Missile is at x = 8. 
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Figure 24. Temperature vs. Longitudinal Index, Displaced One Nodal Point Left of 
Centerline for Round 4 of the F-76 Fire Scenerio. 
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Temperature (K) 


NOTE: First Row: Forward Skin of Missile is at x = 16. Aft Skin of Missile is at x = 19. 





Second Row: Forward Skin of Missile is at x = 9. Aft Skin of Missile is at x = 12. 
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Figure 25. Temperature vs. Longitudinal Index, Displaced One Nodal Point Left of 
Centerline for Round 5 of the F-76 Fire Scenerio. 





V. CONCLUSIONS 


A main objective of this research was to use a commercially-available CFD 
program to develop a three-dimensional model of the thermal effects within a CCL 
compartment. The results could then be compared to previous one-dimensional 
algorithms to achieve a realistic method for non-destructively determining the elapsed time 
for missile cook-off due to a fire in an adjacent compartment. 

The results achieved in these simulations appeared to vary from Callaham’s 
predictions by one order of magnitude. One reason that the one-dimensional model’s 
times are long is that they did not account for reflected radiation received from two of the 
four side bulkheads, the upper and lower steel supports, and off of the other canisters 
within the launcher. On the other hand, the elapsed times to cook-off in this study could 
be short due to the fact that the canisters were modeled as solid titanium cylinders with a 
higher thermal diffusivity than the actual concentric canisters. 

One method to determine the accuracy of the CFD results is by viewing the 
residual plots produced by CFD-ACE during the simulation (Appendix N). The software 
developers at CFDRC recommend a five order of magnitude drop in the residuals of all of 
the calculated variables (i.e., velocity, turbulence quantities, radiation, enthalpy and 
pressure). Appendix N shows that this was only achieved for the turbulence quantities, 
and even then, only in the first time step. The enthalpy managed to drop three orders of 
magnitude, which is not entirely unacceptable. 

The number of input items that could have been changed to achieve this uniform 


drop in residual errors is infinite. For example, changes could have been made in the 
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relaxation values, the size of the time step, the distribution of the nodal points, or even in 
the types of differencing techniques used. However, the decision was made to vary only 
the longitudinal distribution of nodes in the geometric model and leave all other input data 
the same for every simulation. 

Another possible method to determine the accuracy of the data 1s by increasing the 
nodal distribution. This is the motive for the multiple rounds of each single-row and two- 
row fire scenerio. The temperature profiles in Appendices D through M show consistent 
values of temperature distribution and elapsed time to cook-off as the grid distribution is 
increased in the longitudinal direction. 

Although the value for the elapsed time to cook-off, as predicted in this study, may 
not be precise, the figures given in the back of this report clearly show an accurate 
temperature profile for the simulated fire scenerios. The three dimensional profiles 
accurately depict the relative size and direction of the circulation of air generated by the 
hot fire boundary. The planar cutaway and temperature-index plots of these profiles 
provide a visualization of the temperature distribution across the launcher. These pictures 
show the immediate radiative heat transfer occurring on all metal surfaces and the 
relatively slow conductive heat transfer through the titanium canisters. The planar 
cutaway of the temperature distribution also shows how the single-row of canisters is 
cooled by the circulating air flow (i.e., the extended blue area overlapping the canister 
circle), and how this effect 1s negated by the placement of a second row of canisters 
directly behind the first row, with a shift in the convection to the backside of the second 


row. 
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The time-temperature results also depict accurate descriptions of the phenomena 
within the center CCL canister of the first row. Figure 15 shows the increase in the time it 
takes to heat the front of the canister as the second row is added. With only one row in 
the CCL enclosure, the radiation intensity 1s strong on the front skin of the canister wall 
since the radiation does not strike a second surface until it is has reached the forward 
bulkhead of the enclosure. The distance from the aft bulkhead of the launcher to the 
second row of canisters is much closer than it is to the forward bulkhead of the enclosure, 
therefore, the radiation no longer intensely strikes only the front row, but is now able to 
also intensely strike the closer second row. 

Furthermore, Figure 15 shows the decrease in the time it takes to heat the back of 
the center canister of the front row, as a second row is added. With only one row in the 
CCL enclosure, the circulation of air cooled the back skin of the canisters, that were 
receiving reflected radiation from the forward bulkhead of the enclosure. Now, with the 
addition of a second row of canisters, the back skin of the second row receives this 
cooling effect, leaving a relatively stagnant air gap between the two rows. Additionally, 
the reflected radiation off of the front skin of the second row of canisters is much stronger 
on the back skin of the first row, because the distance between the canisters is closer than 
the distance between the back skin of the first row and the forward bulkhead of the 
enclosure (see Figures 8 and 9). 

Appendices O and P contain the results of an additional geometric model of the 
CCL enclosure using the F-76 and Exocet fire scenerios. This model also contained two 


rows of canisters, however, the forward bulkhead of the launcher was placed directly 
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behind the second row of cylinders, at the same distance that separated the aft bulkhead 
from the first row. The purpose of this simulation was to gain an understanding of how 
the temperatures in the two rows of canisters were affected by removing the large volume 
of excess space that the original model contained. 

As expected, Appendices O and P show that the reflected radiation from the 
forward bulkhead of the enclosure increased as the distance between it and the second row 
of canisters was decreased. Additionally, the tighter enclosure reduced the amount of 
circulation within the space, and thus the back skin of the second row of canisters does 
not receive a cooling effect as shown in the original two-row model. The time- 
temperature results plotted for the front row monitored points show that only the front 
skin of the canisters are affected by the decrease in the longitudinal distance between the 
aft and forward bulkheads of the enclosure. The cook-off times for the front skin in the 
reduced geometric launcher are decreased by approximately one-half the values in the full- 
size intermediate-level CCL compartment, due to an increase in the radiation intensity on 
the front skin of the canisters in the first row. 

In general, this study attempted to pioneer the use of today’s CFD technology to 
simulate the effects of shipboard fires, for the development of improved firefighting 
detection and extinguishing systems and procedures. The visual profiles produced by the 
CFD-VIEW program have proven to accurately predict the temperature profiles in the 
modeling of the thermal effects in a CCL caused by a fire in an adjacent compartment. As 
the power of the computer and accuracy of CFD programs improve, this method of 


shipboard fire simulation will one day become the norm for system design and engineering. 
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VI. RECOMMENDATIONS 


The following recommendations are made in the continuation of this study: 


The first step that needs to be taken is to run the software on a faster and more 
powerful computer system, such as a Cray. 


Decrease the time step below 0.37 seconds. Callaham’s algorithm required a time step 
below this value in order to satisfy stability requirements. It would be interesting to 
see how this value effects the three-dimensional CFD model, especially in the residual 
plots. 


Increase the geometry of the model to include the remaining five rows of canisters. 
Add one row at a time to see how each additional row affects the temperature 
distribution and elapsed time to propellant cook-off. 


Completely model a single concentric canister, placing the minimum cook-off 
temperatures on the canister’s skin to determine the effects they have on the internal 
propellant. 


The CFD-ACE package includes a combustion algonthm. Use this algonthm to model 
the effects of an actual shipboard fire, to develop optimum firefighting procedures and 
systems. 


The CFD-ACE package also includes a spray algorithm. Use this algorithm to model 


the effects of sprinkler location, in comparison to the location of possible fire sources 
in a Main engine room or magazine. 
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APPENDIX A. OVERVIEW OF CCL SIMULATIONS 
This matrix shows the five rounds of simulations used to model the Exocet and 
F-76 fire scenerios. The one row model used three increments of nodal distribution in the 


longitudinal direction and the two row model used two increments. 


(A, B) Exocet Fire F-76 Fire 
Round | 
(4, 8)=14 X X 
B | Round 2 
(4, 12) 18 x A 
OOO00O0 Round 3 
o (Ga20) 28 X x 
- Round 4 
B (2,4) =18 A X 


OOOOO|. Round 5 xX x 
OOOOO} ¢3.8)=24 


Figure 26. Overview of CCL Simulations. 
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APPENDIX B. NODAL DISTRIBUTION IN THE 
LONGITUDINAL DIRECTION 





Figure 27. Horizontal Grid Distribution for Round 1. 
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Figure 28. Horizontal Grid Distribution for Round 2. 
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Figure 29. Horizontal Grid Distribution for Round 3. 
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Figure 30. Horizontal Grid Distribution for Round 4. 
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Figure 31. Horizontal Grid Distribution for Round 5. 
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APPENDIX C. CFD-ACE INPUT FILE. 


This is a sample CFD-ACE input file used to run the CCL simulations. All of the 
experiments used the same general format. The only two differences in the input files for 
each of the runs is in the number of longitudinal nodes and the number of canisters (ie. five 
for the single row runs and ten for the two row runs.). This particular file is for an F-76 


fire scenerio involving a single row of canisters. 


MODEL Round _ 3 
* 
GEOMETRY 
GRID 3D BFC 
READ GRID FROM Round_3.PFG 
* Cell Types 
* name: tomahawk 
SOLID titanium 4 6 8 14 20 22 
* name: tomahawk 
SOLID titanium 4 6 1 7 20 22 
* name: tomahawk 
SOLID titanium 10 12 8 14 20 22 
* name: tomahawk 
SOLID titanium 10 12 1 7 20 22 
* name: tomahawk 
SOLID titanium 16 18 8 14 20 22 
* name: tomahawk 
SOLID titanium 16 18 1 7 20 22 
* name: tomahawk 
SOLID titanium 22 24 8 14 20 22 
* name: tomahawk 
SOLID titanium 22 24 1 7 20 22 
* name: tomahawk 
SOLID titanium 28 30 8 14 20 22 
* name: tomahawk 
SOLID titanium 28 30 1 7 20 22 
END 
* 
PROBLEM TYPE 
SOLVE FLOW TURBULENCE HEAT RADIATION 
UNSTEADY TF = 0 TL = 600 STEPS = 60 


ae 


END 
* 
PROPERTIES 
DENSITY INVERSE. T 
VISCOSITY CONSTANT KINEMATIC 1.71e-05 
CONDUCTIVITY PRANDTL 0.7 
SPECIFIC_HEAT 1000 
PRT 1 
SOLID_PROPERTIES titanium K = 21.9 CP = 522 RHO = 4930 EMISS = 0.6 
RADIATION EMISS = 0.82 ABSOR = 0 SCATTER = 0 
END 
€ 
MODELS 
TURBULENCE_MODEL KE 
END 


* 


*** Boundary Conditions *** 

* 

BOUNDARY _ CONDITIONS 

* Gravity boundary condition specified 
GRAV_X0 
GRAV_Y -9.81 
RHOREF AUTO 
GRAV Z0 

* boundary condition: Default 
WALL 11814119 WEST 
U=0V=0W=0 

* boundary condition: top_deck 
WALL 13 14 14119 NORTH 
U=0V=0W=0 

* boundary condition: Default 
WALL 1381411 LOW 
U=0V=0W=0 

* boundary condition: Default 
WALL 1117119 WEST 
U=0V=0W=0 

* boundary condition: Default 
WALL 13 11119 SOUTH 
U=0V=0W=0 

* boundary condition: Default 
WALL 131711 LOW 
U=0V=0W=0 

* boundary condition: top_deck 
WALL 4 6 14 141 19 NORTH 
U=0V=0W=0 
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* boundary condition: Default 
WALL 468 1411 LOW 
U=0V=0W=0 

* boundary condition: Default_ 10 

* boundary condition: Default 
WALL 4611119 SOUTH 
U=0 V=0W=0 

* boundary condition: Default 
WALL 461711 LOW 
U=0V=0W=0 

* boundary condition: Default_ 10 

* boundary condition: top_deck 
WALL 79 14 14119 NORTH 
U=0V=0W=0 

* boundary condition: Default 
WALL 7981411 LOW 
U=0V=0W=0 

* boundary condition: Default 
WALL 7911119 SOUTH 
U=0V=0W=0 

* boundary condition: Default 
WALL 791711 LOW 
U=0V=0W=0 

* boundary condition: top deck 
WALL 10 12 14 14119 NORTH 
U=0V=0W=0 

* boundary condition: Default 
WALL 10 1281411 LOW 
U=0V=0W=0 

* boundary condition: Default_ 10 

* boundary condition: Default 
WALL 10 1211119 SOUTH 
U=0V=0W=0 

* boundary condition: Default 
WALL 10121711 LOW 
U=0V=0W=0 

* boundary condition: Default_10 

* boundary condition: top deck 
WALL 13 15 14 14119 NORTH 
U=0V=0W=0 

* boundary condition: Default 
WALL 13 158 1411 LOW 
U=0V=0W=0 

* boundary condition: Default 
WALL 13 1511119 SOUTH 
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U=0V=0W=0 

* boundary condition: Default 
WALL 13 151711 LOW 
U=0V=0W=0 

* boundary condition: top_deck 
WALL 16 18 14 14119 NORTH 
U=0V=0W=0 

* boundary condition: Default 
WALL 16 188 1411 LOW 
U=0V=0W=0 

* boundary condition: Default_ 10 

* boundary condition: Default 
WALL 16 18 11119 SOUTH 
U=0V=0W=0 

* boundary condition: Default 
WALL 16 181711 LOW 
U=0V=0W=0 

* boundary condition: Default_10 

* boundary condition: top deck 
WALL 19 21 14 14119 NORTH 
U=0V=0W=0 

* boundary condition: Default 
WALL 19 218 1411 LOW 
U=0V=0W=0 

* boundary condition: Default 
WALL 19 21 11119 SOUTH 
U=0V=0W=0 

* boundary condition: Default 
WALL 19211711 LOW 
U=0V=0W=0 

* boundary condition: top deck 
WALL 22 24 1414119 NORTH 
U=0V=0W=0 

* boundary condition: Default 
WALL 22 2481411 LOW 
U=0V=0W=0 

* boundary condition: Default_ 10 

* boundary condition: Default 
WALL 22 2411119 SOUTH 
U=0V=0W=0 

* boundary condition: Default 
WALL 22 241711 LOW 
U=0V=0W=0 

* boundary condition: Default_ 10 

* boundary condition: top deck 
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WALL 25 27 14 14119 NORTH 
U=0V=0W=0 
* boundary condition: Default 
WALL 25 278 1411 LOW 
U=0V=0W=0 
boundary condition: Default 
WALL 25 2711119 SOUTH 
U=0V=0W=0 
boundary condition: Default 
WALL 25 271711 LOW 
U=0V=0W=0 
boundary condition: top_deck 
WALL 28 30 14 141 19 NORTH 
U=0V=0W=0 
boundary condition: Default 
WALL 28 308 1411 LOW 
U=0V=0W=0 
boundary condition: Default_ 10 
boundary condition: Default 
WALL 28 3011119 SOUTH 
U=0V=0W=0 
boundary condition: Default 
WALL 28 301711 LOW 
U=0V=0W=0 
boundary condition: Default_ 10 
boundary condition: Default 
WALL 33 33 8 141 19 EAST 
U=0V=0W=0 
boundary condition: top deck 
WALL 31 33 1414119 NORTH 
U=0V=0W=0 
boundary condition: Default 
WALL 31 33 81411 LOW 
U=0V=0W=0 
boundary condition: Default 
WALL 33 33 17119 EAST 
U=0V=0W=0 
* boundary condition: Default 
WALL 31 33 11119 SOUTH 
U=0V=0W=0 
* boundary condition: Default 
WALL 31 331711 LOW 
U=0V=0W=0 
boundary condition: Default 
WALL 1 18 142022 WEST 


* 


* 


* 


* 


x 


* 


xk 


% 


* 


* 


* 
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b= 0 V=O WwW =0 


* boundary condition: Default_10 


* 


* 


boundary condition: top_ deck 
WALL 1 3 14 14 2022 NORTH 
U=0V=0W=0 

boundary condition: Default 
WALL 1 1 1 72022 WEST 
U=0V=0W=0 


* boundary condition: Default_10 


* 


boundary condition: Default 
WALL 13 112022 SOUTH 
U=0V=0W=0 


* boundary condition: Default_10 
* boundary condition: Default_10 


* 


boundary condition: top_deck 
WALL 4 6 14 1420 22 NORTH 
U=0V=0W=0 


* boundary condition: Wall_ 10 
* boundary condition: Wall 10 


* 


boundary condition: Default 
WALL 461 12022 SOUTH 
U=0V=0W=0 


* boundary condition: Default_10 
* boundary condition: Default_10 


* 


boundary condition: top_deck 
WALL 79 14 142022 NORTH 
U=0V=0W=0 


* boundary condition: Wall_ 10 


* 


* 
* 
* 


* 
* 
* 


* 
* 
ok 


boundary condition: Default 
WALL 79112022 SOUTH 
U=0V=0W=0 

boundary condition: Default_10 
boundary condition: Default_ 10 
boundary condition: top deck 
WALL 10 12 14 1420 22 NORTH 
U=0 V=0W=0 

boundary condition: Wall 10 
boundary condition: Wall 10 
boundary condition: Default 
WALL 10 12112022 SOUTH 
U=0V=0W=0 

boundary condition: Default_10 
boundary condition: Default_ 10 
boundary condition: top_ deck 
WALL 13 15 14 1420 22 NORTH 
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* 
* 


* 
* 
* 


* 
* 
* 


* 
* 
* 


* 
* 


U=0V=0W=0 

boundary condition: Wall_ 10 
boundary condition: Default 
WALL 13 15 1 12022 SOUTH 
U=0V=0W=0 

boundary condition: Default_ 10 
boundary condition: Default_10 
boundary condition: top_deck 
WALL 16 18 14 1420 22 NORTH 
U=0V=0W=0 

boundary condition: Wall_ 10 
boundary condition: Wall_ 10 
boundary condition: Default 
WALL 16 18 1 1 2022 SOUTH 
U=0V=0W=0 

boundary condition: Default_ 10 
boundary condition: Default_ 10 
boundary condition: top_deck 
WALL 19 21 14 1420 22 NORTH 
U=0V=0W=0 

boundary condition: Wall_ 10 
boundary condition: Default 
WALL 19 21 1 12022 SOUTH 
U=0V=0W=0 


* boundary condition: Default_10 


* 
* 


boundary condition: Default_10 
boundary condition: top_ deck 
WALL 22 24 14 14 20 22 NORTH 
U=0V=0W=0 


* boundary condition: Wall_ 10 


* 
* 


* 
* 
* 


* 
* 


* 


* 


boundary condition: Wall 10 
boundary condition: Default 
WALL 22 24 1 12022 SOUTH 
U=0V=0W=0 

boundary condition: Default_ 10 
boundary condition: Default_ 10 
boundary condition: top_deck 
WALL 25 27 14 14 20 22 NORTH 
U=0V=0W=0 

boundary condition: Wall. 10 
boundary condition: Default 
WALL 25 27 1 1 2022 SOUTH 
U=0V=0W=0 

boundary condition: Default_ 10 
boundary condition: Default_ 10 
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* 


ete. woe 


* 


boundary condition: top_deck 
WALL 28 30 14 142022 NORTH 
U=0V=0W=0 

boundary condition: Wall_ 10 
boundary condition: Wall_ 10 
boundary condition: Default 
WALL 28 30 1 1 2022 SOUTH 
U=0V=0W=0 

boundary condition: Default_10 


* boundary condition: Default 


* 


ok 


* 


*% 


* 


* 


* 


* 


% 


*% 


WALL 33 33 8 142022 EAST 
U=0V=0W=0 

boundary condition: top_deck 
WALL 31 33 14 1420 22 NORTH 
U=0V=0W=0 

boundary condition: Default 
WALL 33 33 1 72022 EAST 
U=0V=0W=0 

boundary condition: Default 
WALL 31 33 1 12022 SOUTH 
U=0V=0W=0 

boundary condition: Default 
WALL 1 1 8 1423 27 WEST 
U=0V=0W=0 

boundary condition: top_deck 
WALL 1 3 14 14 23 27 NORTH 
U=0V=0W=0 

boundary condition: Fire_bulkhead 
WALL 13 8 142727 HIGH 
U=0V=0W=0 

boundary condition: Default 
WALL 11172327 WEST 
U=0V=0W=0 

boundary condition: Default 
WALL 13 112327 SOUTH 
U=0V=0W=0 

boundary condition: Fire_bulkhead 
WALL 13 1 72727 HIGH 
U=0V=0W=0T=PROF T 
T 4 

0 100 200 300 

300 655 1010 1365 

boundary condition: top_deck 
WALL 4 6 14 14 23 27 NORTH 
U=0V=0W=0 
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* boundary condition: Fire_bulkhead 
WALL 4 6 8 14 27 27 HIGH 
U=0V=0W=0 

* boundary condition: Default 
WALL 46 1 123 27 SOUTH 
U=0V=0W=0 

* boundary condition: Fire_bulkhead 
WALL 4 6 17 2727 HIGH 
U=0V=0W=0T=PROF T 
T4 
0 100 200 300 
300 655 1010 1365 

* boundary condition: top_deck 
WALL 7 9 14 14 23 27 NORTH 
U=0V=0W=0 

* boundary condition: Fire_bulkhead 
WALL 7 9 8 1427 27 HIGH 
U=0V=0W=0 

* boundary condition: Default 
WALL 79 1 1 23 27 SOUTH 
U=0V=0W=0 

* boundary condition: Fire_bulkhead 
WALL 79 17 2727 HIGH 
U=0V=0W=0T=PROF T 
T 4 
0 100 200 300 
300 655 1010 1365 

* boundary condition: top deck 
WALL 10 12 14 14 23 27 NORTH 
U=0V=0W=0 

* boundary condition: Fire_bulkhead 
WALL 10 12 8 14 2727 HIGH 
U=0V=0W=0 

* boundary condition: Default 
WALL 10 12 1 1 23 27 SOUTH 
U=0V=0W=0 

* boundary condition: Fire bulkhead 
WALL 10 12172727 HIGH 
U=0 V=0W=0T=PROF T 
T4 
0 100 200 300 
300 655 1010 1365 

* boundary condition: top_deck 
WALL 13 15 14 1423 27 NORTH 
U=0V=0W=0 
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* boundary condition: Fire_bulkhead 
WALL 13 15 8 14 2727 HIGH 
U=0V=0W=0 

* boundary condition: Default 
WARE 13 15 112327 SOUTH 
U=0V=0W=0 

* boundary condition: Fire_bulkhead 
WALL 13 15 172727 HIGH 
U=0V=0W=0T=PROF T 
T 4 
0 100 200 300 
300 655 1010 1365 

* boundary condition: top deck 
WALL 16 18 14 14 23 27 NORTH 
U=0V=0W=0 

* boundary condition: Fire_bulkhead 
WALL 16 18 8 142727 HIGH 
U=0V=0W=0 

* boundary condition: Default 
WALL 16 18 1 1 2327 SOUTH 
U=0V=0W=0 

* boundary condition: Fire_bulkhead 
WALL 16 18 1 72727 HIGH 
U=0V=0W=0T=PROF T 
T4 
0 100 200 300 
300 655 1010 1365 

* boundary condition: top deck 
WALL 19 21 14 1423 27 NORTH 
U=0V=0W=0 

* boundary condition: Fire_bulkhead 
WALL 19 21 8 142727 HIGH 
U=0V=0W=0 

* boundary condition: Default 
WALL 19 21 1 12327 SOUTH 
U=0V=0W=0 

* boundary condition: Fire_bulkhead 
WALL 19 21 1 72727 HIGH 
U=0V=0W=0T=PROF T 
T 4 
0 100 200 300 
300 655 1010 1365 

* boundary condition: top deck 
WALL 22 24 14 14 23 27 NORTH 
U=0V=0W=0 
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* boundary condition: Fire_bulkhead 
WALL 22 24 8 14 27 27 HIGH 
U=0V=0W=0 

* boundary condition: Default 
WALL 22 24 1 1 23 27 SOUTH 
U=0V=0W=0 

* boundary condition: Fire_bulkhead 
WALL 22 24 172727 HIGH 
U=0V=0W=0T=PROF T 
T 4 
0 100 200 300 
300 655 1010 1365 

* boundary condition: top_deck 
WALL 25 27 14 14 23 27 NORTH 
U=0V=0W=0 

* boundary condition: Fire bulkhead 
WALL 25 27 8 1427 27 HIGH 
U=0V=0W=0 

* boundary condition: Default 
WALL 25 27 1 123 27 SOUTH 
U=0V=0W=0 

* boundary condition: Fire_bulkhead 
WALL 25 27 172727 HIGH 
U=0 V=0W=0 T= PROF T 
T 4 
0 100 200 300 
300 655 1010 1365 

* boundary condition: top_deck 
WALL 28 30 14 14 23 27 NORTH 
U=0V=0W=0 

* boundary condition: Fire_bulkhead 
WALL 28 30 8 14 27 27 HIGH 
U=0V=0W=0 

* boundary condition: Default 
WALL 28 30 1 1 23 27 SOUTH 
U=0V=0W=0 

* boundary condition: Fire_bulkhead 
WALL 28 30 1 72727 HIGH 
U=0 V=0W=0T=PROF T 
T4 
0 100 200 300 
300 655 1010 1365 

* boundary condition: Default 
WALL 33 33 8 14 23 27 EAST 
U=0V=0W=0 
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* boundary condition: top_deck 
WALL 31 33 14 1423 27 NORTH 
U=0V=0W=0 

* boundary condition: Fire_bulkhead 
WALL 31 33 8 1427 27 HIGH 
U=0V=0W=0 

* boundary condition: Default 
WALL 33 33 1723 27 EAST 
U=0V=0W=0 

* boundary condition: Default 
WALL 31 33 1 123 27 SOUTH 
U=0V=0W=0 

* boundary condition: Fire_bulkhead 
WALL 31 33 172727 HIGH 
U=0 V=0W=0T=PROF T 
T 4 
0 100 200 300 
300 655 1010 1365 

END 

* 

INITIAL _ CONDITIONS 

* Full field initial conditions 
U=0 V=0W=0 P=0T=300 K =5 D=25L=0 

END 

* 

2K OK OK OK 2 OK OK OK OK OK KOK KK KK KOK KKK Solve 2K OK OK OK OK OK OK 2K OK OK Ko OK OK KOK 2 Ok OK OK OK KK 

* 

*k 


SOLUTION CONTROL 
ALGORITHM SIMPLEC 
S SCHEME UPWIND U V W RHO K DH 
T_SCHEME EULER 
ITERATIONS 100 
C_ITERATIONS 1 
SOLVER WHOLE_JU VW PPKD 
SOLVER CG H 
S ITERATIONS 8 U V W 
S ITERATIONS 30 PP 
S ITERATIONS 5 K D 
S ITERATIONS 1000 H 
INERTIAL FACTOR 1 UVW 
INERTIAL FACTOR 0.2 K D 
INERTIAL FACTOR 0.1 H 
RELAX 0.8 P 
RELAX | RHO 
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RELAX 0.3 T 
RELAX 1 VIS 
RELAX 0.1 RAD 
MINVAL -le+20 U V W 
MINVAL -let+20 P 
MINVAL le-06 RHO 
MINVAL le-10 T VIS 
MINVAL le-10 K D 
MINVAL -le+20 H 
MAXVAL 1le+20 U V W 
MAXVAL 1le+20 P RHO 
MAXVAL 5000 T 
MAXVAL 1le+20 VIS 
MAXVAL 1le+20 K D 
MAXVAL 1le+20 H 
END 
* 
OUTPUT 
PLOT3D ON FORMATTED 
SCALAR FILE 1RHOP TKD 
DIAGNOSTICS OFF 
* inside left of far tube wall 
MONITOR 1 173 21 T 
* inside right of far tube wall 
MONITOR 1 183 21 T 
* inside left of near tube wall 
MONITOR 1 17 3 22 T 
* inside right of near tube wall 
MONITOR 1 18 3 22 T 
* left of far tube wall 
MONITOR 1 17 3 20 T 
* right of far tube wall 
MONITOR 1 18 3 20 T 
* left of near tube wall 
MONITOR 1 17 3 23 T 
* right of near tube wall 
MONITOR | 18 3 23 T 
TIME SAVE 6 
UNIQUE_NAME ON 
END 
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APPENDIX D. EXOCET FIRE - ROUND 1 
The following data corresponds to the round 1 simulation of the Exocet fire 
scenerio, including: 


(1) 3-dimensional temperature profile at a time of 1 minute. 
(2) 3-dimensional temperature profile at a time of 10 minutes. 
(3) Horizontal temperature profile at a time of 1 minute. 

(4) Horizontal temperature profile at a time of 10 minutes. 
(5) Temperature vs. index profile at a time of 1 minute. 

(6) Temperature vs. index profile at a time of 10 minutes. 

(7) Temperature vs. time plot for monitor point 2. 

(8) Temperature vs. time plot for monitor point 3. 
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Figure 32. Exocet Fire - Round 1. 3-Dimensional Temperature Profile 


at a Time of 1 Minute. 
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Figure 33. Exocet Fire - Round 1. 3-Dimensional Temperature Profile 


at a Time of 10 Minutes. 
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Figure 34. Exocet Fire - Round 1. Horizontal Temperature Profile 
at a Time of 1 Minute. 
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Figure 35. Exocet Fire - Round 1. Horizontal Temperature Profile 
at a Time of 10 Minutes. 


77 





Temperature (K) 


1900 
1800 
1790 
1600 
1500 
1400 
1300 
12900 
1100 
1000 
300 
800 
700 
6900 
500 
490 
300 








QO 2 


. 
ee ere eA em a 

} . . 

’ ; 

; j i 5 ’ 

i r i { t 3 : , 

: ! : i : , t : 

H } , t 4 : 5 
SE a a a a th mm fr ee = ny er ert ee re he ee epee Se Se eee eee + 

t ‘ } ' r 

: ' $ ; ; { , ’ : : 

; 
f ; 1 t } ; . 


: , 
: ’ , ‘ ’ ¢ 
rn et ret arene rman me are a raat er a a tt a ae a em ee ee nee ee eet en ere eee ee 
) 3 > 5 , 
’ : : Fe t : 


‘ 
¢ i : i ’ : 1 uy ; ; $ 
» 5 d } ’ 1 ' ; t , b 

‘ . ’ , 

1 ' i ‘ { : ’ , 
° $ 
s ’ 
q 


¢ 
' 


ed 


\ | ; 
; { i 
. 3s 


[ t ; 


are 


: < . 
’ ‘ : , , : 
' ’ : ' * > i 
7 4 } { ‘ * ; ; 
nei ee er ee ere rn en mar er ere pe a ee fi | RS rs © Ate «me 
, ae} , 1 ; 
’ > 5 , ' ‘ 
7 , ’ 








; , 
ee a ede 
: ‘ } 


moe we 


: % 
‘ ; { ’ i ‘ | 
i i f i ? { ‘ , 
ar err nr en ner ae nme monomer te 
| 
2 ; 4 3 ¥ ! H , i ’ 
. i ; , : 
Renee Mess ob 2 sa ee ee eee, 
: j : a i : ‘ ' ; : } 
: ‘ t 1 i ‘ { { 
} 4 : { 
ne at ee em se te a eps I a 8 Seer eREROT pore ate ys Smee aoa oe a ree ee wane = be a os erate plea = aaa Sry eono ~4 wa fee er ee 
1 : : ’ £ 
; : : ' : 
3 t i ; } 
a Fee ee te ere teed ae a ara rare ee ee ee ee ee oe ee ee SS aeateaeneneentl 
§ i ? ‘ t \ 
$ ! : y ’ ; ‘ i 
t ; ; 3 ‘ 
; al 
: 7 a, Ronee ort _=_ =—-. 
; ’ i 1 : ; j : 
' : - ' ' $ : t 
Sn ee te em ee fa Vm: en eee : ey 
4 j , f t t : ¢ 
‘ ' ‘ : ; , 
: : i ; 


, : : 
. e ! : , y ’ 
cea ae ae ee ame a ae Tne MS Ne en gaeeemn eet meee ye oe ee pe ee i el 
- . 3 i ’ s } 
. ’ , 


g ’ 
; > , i ) : 
S } ! j ' } 1 £ 
rere rr a pe a ttre 2 ee eee ee ee erat 
; ' ? : s . 
; , tr : t 
i i i g if $ : 
oo watd>! ee! fe ans ee fa i Le 00 Oe! em am EY © By. ove amen - ei we Bee Se a “he es pape ea ——— maP eas 
. . 2 * ’ ’ 
; : ; ? : } : : ’ 4 
5 ‘ ; é F ‘ c ’ 
{ ‘ : ‘ f 
a oe eee. ee te, ene te ema ne cme mnne me, Se ne en oy: eee i et ee re mm em on oe em meee ee ee mare me meee one a ~ me wee 2 see my 
; ; : t ; : : : 
I , ’ , . 
> ' ‘ t 
t : é i 4 
a c 


U4 
p 
in 
m™ 
~J 
Oo 
O 
O 


17 12 


Longitudinal Index Left of the Centerline on Center Canister. 


Figure 36. Exocet Fire - Round 1. Temperature vs. Index Profile 
at a Time of 1 Minute. 
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Longitudinal Index Left of the Centerline on Center Canister. 


Figure 37. Exocet Fire - Round 1. Temperature vs. Index Profile 
at a Time of 10 Minutes. 
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Figure 38. Exocet Fire - Round 1. Temperature vs. Time Plot for Monitor Point 2. 
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Figure 39. Exocet Fire - Round 1. Temperature vs. Time Plot for Monitor Point 3. 
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APPENDIX E. EXOCET FIRE - ROUND 2 
The following data corresponds to the round 2 simulation of the Exocet fire 
scenerio, including: 


(1) 3-dimensional temperature profile at a time of one minute. 
(2) 3-dimensional temperature profile at a time of 10 minutes. 
(3) Horizontal temperature profile at a time of one minute. 
(4) Horizontal temperature profile at a time of 10 minutes. 
(5) Temperature vs. index profile at a time of one minute. 

(6) Temperature vs. index profile at a time of 10 minutes. 

(7) Temperature vs. time plot for monitor point 1. 

(8) Temperature vs. time plot for monitor point 2. 

(9) Temperature vs. time plot for monitor point 3. 
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Figure 40. Exocet Fire - Round 2. 3-Dimensional Temperature Profile 


at a Time of 1 Minute. 
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Figure 41. Exocet Fire - Round 2. 3-Dimensional Temperature Profile 
at a Time of 10 Minutes. 
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Figure 42. Exocet Fire - Round 2. Horizontal Temperature Profile 
at a Time of 1 Minute. 
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Figure 43. Exocet Fire - Round 2. Horizontal Temperature Profile 
at a Time of 10 Minutes. 
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Longitudinal Index Left of the Centerline on Center Canister. 
Figure 44. Exocet Fire - Round 2. Temperature vs. Index Profile 
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Longitudinal Index Left of the Centerline on Center Canister. 
Figure 45. Exocet Fire - Round 2. Temperature vs. Index Profile 
at a Time of 10 Minutes. 
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Figure 46. Exocet Fire - Round 2. Temperature vs. Time Plot for Monitor Point 1. 
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Figure 47. Exocet Fire - Round 2. Temperature vs. Time Plot for Monitor Point 2. 
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Figure 48. Exocet Fire - Round 2. Temperature vs. Time Plot for Monitor Point 3. 
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APPENDIX F. EXOCET FIRE - ROUND 3 
The following data corresponds to the round 3 simulation of the Exocet fire 
scenerio, including: 


(1) 3-dimensional temperature profile at a time of one minute. 
(2) 3-dimensional temperature profile at a time of 10 minutes. 
(3) Horizontal temperature profile at a time of one minute. 
(4) Horizontal temperature profile at a time of 10 minutes. 
(5) Temperature vs. index profile at a time of one minute. 

(6) Temperature vs. index profile at a time of 10 minutes. 

(7) Temperature vs. time plot for monitor point 1. 

(8) Temperature vs. time plot for monitor point 2. 

(9) Temperature vs. time plot for monitor point 3. 
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Figure 49. Exocet Fire - Round 3. 3-Dimensional Temperature Profile 


at a Time of 1 Minute. 
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Figure 50. Exocet Fire - Round 3. 3-Dimensional Temperature Profile 
at a Time of 10 Minutes. 
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Figure 51. Exocet Fire - Round 3. Horizontal Temperature Profile 
at a Time of 1 Minute. 
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Figure 52. Exocet Fire - Round 3. Horizontal Temperature Profile 


at a Time of 10 Minutes. 
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Figure 53. Exocet Fire - Round 3. Temperature vs. Index Profile 
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Figure 54. Exocet Fire - Round 3. Temperature vs. Index Profile 
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Figure 55. Exocet Fire - Round 3. Temperature vs. Time Plot for Monitor Point 1. 
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Figure 56. Exocet Fire - Round 3. Temperature vs. Time Plot for Monitor Point 2. 
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Figure 57. Exocet Fire - Round 3. Temperature vs. Time Plot for Monitor Point 3. 


102 


APPENDIX G. EXOCET FIRE - ROUND 4 
The following data corresponds to the round 4 simulation of the Exocet fire 
scenerio, including: 


(1) 3-dimensional temperature profile at a time of one minute. 
(2) 3-dimensional temperature profile at a time of 10 minutes. 
(3) Horizontal temperature profile at a time of one minute. 
(4) Horizontal temperature profile at a time of 10 minutes. 
(5) Temperature vs. index profile at a time of one minute. 

(6) Temperature vs. index profile at a time of 10 minutes. 

(7) Temperature vs. time plot for monitor point 1. 

(8) Temperature vs. time plot for monitor point 2. 

(9) Temperature vs. time plot for monitor point 3. 
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Figure 58. Exocet Fire - Round 4. 3-Dimensional Temperature Profile 
at a Time of 1 Minute. 
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Figure 59. Exocet Fire - Round 4. 3-Dimensional Temperature Profile 
at a Time of 10 Minutes. 
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Figure 60. Exocet Fire - Round 4. Horizontal Temperature Profile 
at a Time of | Minute. 


106 





1936.0000 
1767.0000 
1604.0000} 
1441.00004 
1278.0000 





re 


1115.00004 
957.0000 fe 


ad 


+ 
o 
oS 


789.0000 4 
626.0000- 
463.0000- 
590.8000 





Figure 61. Exocet Fire - Round 4. Horizontal Temperature Profile 
at a Time of 10 Minutes. 
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Figure 62. Exocet Fire - Round 4. Temperature vs. Index Profile 


at a Time of 1 Minute. 
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Figure 63. Exocet Fire - Round 4. Temperature vs. Index Profile 


at a Time of 10 Minutes. 
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Figure 64. Exocet Fire - Round 4. Temperature vs. Time Plot for Monitor Point 1. 
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Figure 65. Exocet Fire - Round 4. Temperature vs. Time Plot for Monitor Point 2. 
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Figure 66. Exocet Fire - Round 4. Temperature vs. Time Plot for Monitor Point 3. 
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APPENDIX H. EXOCET FIRE - ROUND 5 
The following data corresponds to the round 5 simulation of the Exocet fire 
scenerio, including: 


(1) 3-dimensional temperature profile at a time of one minute. 
(2) 3-dimensional temperature profile at a time of 10 minutes. 
(3) Horizontal temperature profile at a time of one minute. 
(4) Horizontal temperature profile at a time of 10 minutes. 
(5) Temperature vs. index profile at a time of one minute. 

(6) Temperature vs. index profile at a time of 10 minutes. 

(7) Temperature vs. time plot for monitor point 1. 

(8) Temperature vs. time plot for monitor point 2. 

(9) Temperature vs. time plot for monitor point 3. 
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Figure 67. Exocet Fire - Round 5. 3-Dimensional Temperature Profile 
at a Time of 1 Minute. 
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Figure 68. Exocet Fire - Round 5. 3-Dimensional Temperature Profile 
at a Time of 10 Minutes. 
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Figure 69. Exocet Fire - Round 5. Horizontal Temperature Profile 


at a Time of 1 Minute. 
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Figure 70. Exocet Fire - Round 5. Horizontal Temperature Profile 


at a Time of 10 Minutes. 
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Figure 71. Exocet Fire - Round 5. Temperature vs. Index Profile 


at a Time of {1 Minute. 
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Figure 72. Exocet Fire - Round 5. Temperature vs. Index Profile 


at a Time of 10 Minutes. 
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Figure 73. Exocet Fire - Round 5. Temperature vs. Time Plot for Monitor Point 1. 
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Figure 74. Exocet Fire - Round 5. Temperature vs. Time Plot for Monitor Point 2. 
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Figure 75. Exocet Fire - Round 5. Temperature vs. Time Plot for Monitor Point 3. 


APPENDIX I. F-76 FIRE - ROUND 1 
The following data corresponds to the round 1 simulation of the F-76 fire scenerio, 
including: 


(1) 3-dimensional temperature profile at a time of 1.67 minutes. 
(2) 3-dimensional temperature profile at a time of 10 minutes. 
(3) Horizontal temperature profile at a time of 1.67 minutes. 
(4) Horizontal temperature profile at a time of 10 minutes. 

(5) Temperature vs. index profile at a time of 1.67 minutes. 

(6) Temperature vs. index profile at a time of 10 minutes. 

(7) Temperature vs. time plot for monitor point 1. 

(8) Temperature vs. time plot for monitor point 2. 

(9) Temperature vs. time plot for monitor point 3. 
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Figure 76. F-76 Fire - Round 1. 3-Dimensional Temperature Profile 
at a Time of 1.67 Minutes. 
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Figure 77. F-76 Fire - Round 1. 3-Dimensional Temperature Profile 
at a Time of 10 Minutes. 
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Figure 78. F-76 Fire - Round 1. Horizontal Temperature Profile 
at a Time of 1.67 Minutes. 
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Figure 79. F-76 Fire - Round 1. Horizontal Temperature Profile 


at a Time of 10 Minutes. 
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Figure 80. F-76 Fire - Round 1. Temperature vs. Index Profile 


at a Time of 1.67 Minutes. 
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Figure 81. F-76 Fire - Round 1. Temperature vs. Index Profile 


ata Time of 10 Minutes. 
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Figure 82. F-76 Fire - Round 1. Temperature vs. Time Plot for Monitor Point 1. 
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Figure 83. F-76 Fire - Round 1. Temperature vs. Time Plot for Monitor Point 2. 
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Figure 84. F-76 Fire - Round 1. Temperature vs. Time Plot for Monitor Point 3. 


APPENDIX J. F-76 FIRE - ROUND 2 
The following data corresponds to the round 2 simulation of the F-76 fire scenerio, 
including: 


(1) 3-dimensional temperature profile at a time of one minute. 
(2) 3-dimensional temperature profile at a time of 10 minutes. 
(3) Horizontal temperature profile at a time of one minute. 
(4) Horizontal temperature profile at a time of 10 minutes. 
(5) Temperature vs. index profile at a time of one minute. 

(6) Temperature vs. index profile at a time of 10 minutes. 

(7) Temperature vs. time plot for monitor point 1. 

(8) Temperature vs. time plot for monitor point 2. 

(9) Temperature vs. time plot for monitor point 3. 
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Figure 85. F-76 Fire - Round 2. 3-Dimensional Temperature Profile 
at a Time of I Minute. 
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Figure 86. F-76 Fire - Round 2. 3-Dimensional Temperature Profile 
at a Time of 10 Minutes. 
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Figure 87. F-76 Fire - Round 2. Horizontal Temperature Profile 
at a Time of 1 Minute. 
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Figure 88. F-76 Fire - Round 2. Horizontal Temperature Profile 
at a Time of 10 Minutes. 
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Figure 89. F-76 Fire - Round 2. Temperature vs. Index Profile 


at a Time of 1 Minute. 
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Figure 90. F-76 Fire - Round 2. Temperature vs. Index Profile 


at a Time of 10 Minutes. 
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Figure 91. F-76 Fire - Round 2. Temperature vs. Time Plot for Monitor Point 1. 
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Figure 92. F-76 Fire - Round 2. Temperature vs. Time Plot for Monitor Point 2. 
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Figure 93. F-76 Fire - Round 2. Temperature vs. Time Plot for Monitor Point 3. 


APPENDIX K. F-76 FIRE - ROUND 3 
The following data corresponds to the round 3 simulation of the F-76 fire scenerio, 
including: 


(1) 3-dimensional temperature profile at a time of one minute. 
(2) 3-dimensional temperature profile at a time of 10 minutes. 
(3) Horizontal temperature profile at a time of one minute. 
(4) Horizontal temperature profile at a time of 10 minutes. 
(5) Temperature vs. index profile at a time of one minute. 

(6) Temperature vs. index profile at a time of 10 minutes. 

(7) Temperature vs. time plot for monitor point 1. 

(8) Temperature vs. time plot for monitor point 2. 

(9) Temperature vs. time plot for monitor point 3. 
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Figure 94. F-76 Fire - Round 3. 3-Dimensional Temperature Profile 
at a Time of 1 Minute. 


144 





1370.0000 
1263.0600 
1156.0600 
1049.0000 
942.0000 
835.0000 
728.0000 
621.0000 
514.0000 
407.0000 
500.0600 





Figure 95. F-76 Fire - Round 3. 3-Dimensional Temperature Profile 
at a Time of 10 Minutes. 
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Figure 96. F-76 Fire - Round 3. Horizontal Temperature Profile 


at a Time of 1 Minute. 
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Figure 97, F-76 Fire - Round 3. Horizontal Temperature Profile 
at a Time of 10 Minutes. 
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Longitudinal Index Left of the Centerline on Center Canister. 


Figure 98. F-76 Fire - Round 3. Temperature vs. Index Profile 
at a Time of I Minute. 
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Figure 99. F-76 Fire - Round 3. Temperature vs. Index Profile 


at a Time of 10 Minutes. 
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Figure 100. F-76 Fire - Round 3. Temperature vs. Time Plot for Monitor Point 1. 
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Figure 101. F-76 Fire - Round 3. Temperature vs. Time Plot for Monitor Point 2. 
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Figure 102. F-76 Fire - Round 3. Temperature vs. Time Plot for Monitor Point 3. 


APPENDIX L. F-76 FIRE - ROUND 4 
The following data corresponds to the round 4 simulation of the F-76 fire scenerio, 
including: 


(1) 3-dimensional temperature profile at a time of one minute. 
(2) 3-dimensional temperature profile at a time of 10 minutes. 
(3) Horizontal temperature profile at a time of one minute. 
(4) Horizontal temperature profile at a time of 10 minutes. 
(5) Temperature vs. index profile at a time of one minute. 

(6) Temperature vs. index profile at a time of 10 minutes. 

(7) Temperature vs. time plot for monitor point 1. 

(8) Temperature vs. time plot for monitor point 2. 

(9) Temperature vs. time plot for monitor point 3. 
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Figure 103. F-76 Fire - Round 4. 3-Dimensional Temperature Profile 
at a Time of 1 Minute. 
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Figure 104. F-76 Fire - Round 4. 3-Dimensional Temperature Profile 
at a Time of 10 Minutes. 
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Figure 105, F-76 Fire - Round 4. Horizontal Temperature Profile 
at a Time of 1 Minute. 
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Figure 106. F-76 Fire - Round 4. Horizontal Temperature Profile 
at a Time of 10 Minutes. 
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Figure 107. F-76 Fire - Round 4. Temperature vs. Index Profile 


at a Time of 1 Minute. 
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Figure 108. F-76 Fire - Round 4. Temperature vs. Index Profile 
at a Time of 10 Minutes. 
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Figure 109. F-76 Fire - Round 4. Temperature vs. Time Plot for Monitor Point 1. 
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Figure 110. F-76 Fire - Round 4. Temperature vs. Time Plot for Monitor Point 2. 
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Figure 111. F-76 Fire - Round 4. Temperature vs. Time Plot for Monitor Point 3. 
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APPENDIX M. F-76 FIRE - ROUND 5 
The following data corresponds to the round 5 simulation of the F-76 fire scenerio, 
including: 


(1) 3-dimensional temperature profile at a time of one minute. 
(2) 3-dimensional temperature profile at a time of 10 minutes. 
(3) Horizontal temperature profile at a time of one minute. 
(4) Horizontal temperature profile at a time of 10 minutes. 
(5) Temperature vs. index profile at a time of one minute. 

(6) Temperature vs. index profile at a time of 10 minutes. 

(7) Temperature vs. time plot for monitor point 1. 

(8) Temperature vs. time plot for monitor point 2. 

(9) Temperature vs. time plot for monitor point 3. 








Figure 112. F-76 Fire - Round 5. 3-Dimensional Temperature Profile 
at a Time of 1 Minute. 
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Figure 113. F-76 Fire - Round 5. 3-Dimensional Temperature Profile 
at a Time of 10 Minutes. 
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Figure 114. F-76 Fire - Round 5. Horizontal Temperature Profile 
at a Time of | Minute. 
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Figure 115. F-76 Fire - Round 5. Horizontal Temperature Profile 
at a Time of 10 Minutes. 
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Figure 116. F-76 Fire - Round 5. Temperature vs. Index Profile 


at a Time of 1 Minute. 
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Figure 117. F-76 Fire - Round 5. Temperature vs. Index Profile 


at a Time of 10 Minutes. 
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Figure 118. F-76 Fire - Round 5. Temperature vs. Time Plot for Monitor Point 1. 
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Figure 119. F-76 Fire - Round 5. Temperature vs. Time Plot for Monitor Point 2. 
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Figure 120. F-76 Fire - Round 5. Temperature vs. Time Plot for Monitor Point 3. 
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APPENDIX N. RESIDUAL ERROR PROFILES. 


The following residual error figures are representative of all Exocet and F-76 fire 
scenerios at time steps of | (i.e., t = 10 sec) and 60 (1.e., t= 10 min). The variables are: 

u = velocity component in the I-direction. 

v = velocity component in the J-direction. 

w = velocity component in the K-direction. 

PP = pressure. 

K = turbulent kinetic energy. 

D = dissipation rate of turbulent kinetic energy. 

H = enthalpy. 


RX = radiation. 


3 





Figure 121. Residual Profile at Time Step 1 (10 seconds). 





Figure 122. Residual Profile at Time Step 60 (10 minutes). 
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APPENDIX O. EXOCET FIRE - CLOSE WALL 
The following data corresponds to a close wall simulation of the Exocet fire 
scenerio, including: 


(1) Horizontal Grid Distribution. 

(2) 3-dimensional temperature profile at a time of one minute. 
(3) 3-dimensional temperature profile at a time of 10 minutes. 
(4) Horizontal temperature profile at a time of one minute. 
(5) Horizontal temperature profile at a time of 10 minutes. 
(6) Temperature vs. index profile at a time of one minute. 

(7) Temperature vs. index profile at a time of 10 minutes. 

(8) Temperature vs. time plot for monitor point 1. 

(9) Temperature vs. time plot for monitor point 2. 

(10) Temperature vs. time plot for monitor point 3. 
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Figure 123. Horizontal Grid Distribution for Close Wall. 
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Figure 124. Exocet Fire - Close Wall. 3-Dimensional Temperature Profile 
at a Time of 1 Minute. 
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Figure 125. Exocet Fire - Close Wall. 3-Dimensional Temperature Profile 
at a Time of 10 Minutes. 
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Figure 126. Exocet Fire - Close Wall. Horizontal Temperature Profile 


at a Time of 1 Minute. 
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Figure 127. Exocet Fire - Close Wall. Horizontal Temperature Profile 


at a Time of 10 Minutes. 
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Longitudinal Index Left of Centerline of Center Canister 


Figure 128. Exocet Fire - Close Wall. Temperature vs. Index Profile 
ata Time of 1 Minute. 
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Longitudinal Index Left of Centerline of Center Canister 


Figure 129. Exocet Fire - Close Wall. Temperature vs. Index Profile 
at a Time of 10 Minutes. 
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Figure 130. Exocet Fire - Close Wall. Temperature vs. Time Plot 
for Monitor Point lL. 
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Figure 131. Exocet Fire - Close Wall. Temperature vs. Time Plot 
for Monitor Point 2. 
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Figure 132. Exocet Fire - Close Wall. Temperature vs. Time 
Plot for Monitor Point 3. 
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APPENDIX P. F-76 FIRE - CLOSE WALL 
The following data corresponds to a close wall simulation of the F-76 fire scenerio, 
including: 


(1) Horizontal Grid Distribution. 

(2) 3-dimensional temperature profile at a time of one minute. 
(3) 3-dimensional temperature profile at a time of 10 minutes. 
(4) Horizontal temperature profile at a time of one minute. 
(5) Horizontal temperature profile at a time of 10 minutes. 
(6) Temperature vs. index profile at a time of one minute. 

(7) Temperature vs. index profile at a time of 10 minutes. 

(8) Temperature vs. time plot for monitor point 1. 

(9) Temperature vs. time plot for monitor point 2. 

(10) Temperature vs. time plot for monitor point 3. 
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Figure 133. Horizontal Grid Distribution for Close Wall. 
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Figure 134. F-76 Fire - Close Wall. 3-Dimensional Temperature Profile 
at a Time of 1 Minute. 
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Figure 135. F-76 Fire - Close Wall. 3-Dimensional Temperature Profile 


at a Time of 10 Minutes. 
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Figure 136. F-76 Fire - Close Wall. Horizontal Temperature Profile 
at a Time of 1 Minute. 


[9s 


513.0000 
491.6000 
470.2000- 
44B.8000 
427.4000 
406.0000 















384.6000 
363.2000 
341.8000 
320.4000 
299.0000 





T 
1370.0000 


1263.0000 

1156.0000 

1049.0000 
942.0000 
835.0000- 
728.0000 
621.0000 
514.0000 
407.0000 
300.0000 





Figure 137. F-76 Fire - Close Wall. Horizontal Temperature Profile 
at a Time of 10 Minutes. 
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Longitudinal Index Left of Centerline of Center Canister 


Figure 138. F-76 Fire - Close Wall. Temperature vs. Index Profile 
ata Time of 1 Minute. 
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Figure 139. F-76 Fire - Close Wall. Temperature vs. Index Profile 
at a Time of 10 Minutes. 
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Figure 140. F-76 Fire - Close Wall. Temperature vs. Time Plot for Monitor Point 1. 
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Figure 141. F-76 Fire - Close Wall. Temperature vs. Time Plot for Monitor Point 2. 
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Figure 142. F-76 Fire - Close Wall. Temperature vs. Time Plot for Monitor Point 3. 
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